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ABSTRACT 

The drive towards high-work turbines has led to 
designs which can be compact, transonic, supersonic, 
counter rotating, or use a dense drive gas. These 
aggressive designs can lead to strong unsteady 
secondary flows and flow separation. The amplitude 
and extent of these unsteady flow phenomena can be 
amplified at off-design operating conditions. Pre-test 
off-design predictions have been performed for a new 
two-stage supersonic turbine design that is currently 
being tested in air. The simulations were performed 
using a three-dimensional unsteady Navier-Stokes 
analysis, and the predicted results have been 
compared with solutions from a validated meanline 
analysis. 

NOMENCLATURE 

Co - spouting velocity (ft/sec) 

Fr - radial force (Ibf) 


Ft - tangential (circumferential) force (lbf) 

Fx - axial force (lbf) 
m - mass flow (lbm/sec) 

M(abs) - absolute Mach number 
M(rel) - relative Mach number 
P - static pressure (psia) 

Power - Aerodynamic horsepower (HP) 

PR - turbine pressure ratio (total-to-static) 
Pt(abs) - absolute total pressure (psia) 

Pt(rel) - relative total pressure (psia) 

R - reaction based on pressure 

Tt(abs) - absolute total temperature (R) 

Tt(rel) - relative total temperature (R) 

U - Pitchline rotor velocity (ft/sec) 

U/Co - Isentropic velocity ratio 

W - work based on mass-avg (BTU/lbm) 

a - absolute circumferential flow angle, (deg) 

(3 - relative circumferential flow angle, (deg) 

r|(ts) - stage total-to-static efficiency 

T|(tt) - stage total-to-total efficiency 

T](ots) - overall turbine total-to-static efficiency 
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r|(ott) - overall turbine total-to-total efficiency 

to - RPM 


INTRODUCTION 

Modern high-work turbines can be compact, 
transonic, supersonic, counter-rotating, and can use a 
dense drive gas. The vast majority of modem rocket 
turbine designs fall into these categories. These 
turbines are often characterized by large amounts of 
flow unsteadiness. The flow unsteadiness can have a 
major impact on turbine performance and durability. 
For example, the Space Transportation Main Engine 
(STME) fuel turbine, a high-work transonic design 
was found to have an unsteady inter-row shock that 
reduced efficiency by 2 points and increased dynamic 
loading by 24 percent. The Revolutionary Reusable 
Technology Turbopump (RRTT), which uses full 
flow oxygen for its drive gas, was found to shed 
vortices with such energy as to raise serious blade 
durability concerns. In both cases, the sources of the 
problems were uncovered before turbopump testing 
with the application of validated unsteady 
computational fluid dynamics (CFD) to the designs. 
In the case of the RRTT and the Alternate 
Turbopump Development (ATD) turbines, the 
unsteady CFD codes were used not just to identify 
problems, but also to guide designs that mitigate 
problems due to unsteadiness. The incorporation of 
unsteady flow analyses into the design process has led 
to turbine designs with higher performance and fewer 
dynamics problems. References 1-4 are examples of 
the application of unsteady CFD to rocket turbine 
designs. 

Recently, CFD has been used in the design of a two- 
stage supersonic turbine that is currently being tested 
at NASA Marshall Space Flight Center (see Refs. 5 
and 6). Meanline, two-dimensional CFD and three- 
dimensional CFD analyses were used in conjunction 
with optimization techniques to design both the flow 
path and the airfoil geometries. The turbine is being 
tested at both design and off-design operating 
conditions. As part of the pre-test effort unsteady 
three-dimensional flow simulations have been 
performed at several off-design conditions and the 
results compared with results from a validated 
meanline analysis. The objective of the current work 
is to quantify the performance of the turbine at off- 
design flow conditions, as well as to characterize 
changes in the unsteadiness as a function of flow 
condition. 


NUMERICAL PROCEDURES 

The meanline code uses a combination of the one- 
dimensional equations of motion and empirical loss 
models to predict the flow and performance quantities 
in the turbine. The meanline analysis used in this 
study has been anchored on many experimental data 
sets (e.g., subsonic, transonic and supersonic), 
including those shown in Fig. 1. 

The governing equations in the CFD analysis are the 
time-dependent, three-dimensional Reynolds- 
averaged Navier-Stokes equations. The algorithm 
consists of a time marching, implicit, finite-difference 
scheme. The procedure is third-order spatially 
accurate and second-order temporally accurate. The 
inviscid fluxes are discretized according to the 
scheme developed by Roe [7]. The viscous fluxes are 
calculated using standard central differences. An 
approximate-factorization technique is used to 
compute the time rate changes in the primary 
variables. Newton sub-iterations are used at each 
global time step to increase stability and reduce 
linearization errors. For all cases investigated in this 
study, one Newton sub-iteration was performed at 
each time step. The turbulent viscosity is calculated 
using the two-layer Baldwin-Lomax algebraic 
turbulence model [8]. The code has been parallelized 
using Message Passing Interface (MPI) and OpenMP 
application program interfaces (API’s) to reduce the 
computation time for large-scale three-dimensional 
simulations. 

The Navier-Stokes analysis uses O- and H-type zonal 
grids to discretize the flow field and facilitate relative 
motion of the rotating components (see Fig. 2). The 
O-grids are body-fitted to the surfaces of the airfoils 
and generated using an elliptic equation solution 
procedure. They are used to properly resolve the 
viscous flow in the blade passages and to easily apply 
the algebraic turbulence model. The algebraically 
generated H-grids are used to discretize the remainder 
of the flow field. 

The CFD analysis has been validated on several 
supersonic turbine geometries (e.g.. Refs. 1-3). 

BOUNDARY CONDITIONS 

The theory of characteristics is used to determine the 
boundary conditions at the inlet and exit of the 
computational domain. For subsonic inlet flow the 
total pressure, total temperature, and the 
circumferential and radial flow angles are specified as 
a function of the radius. The upstream running 
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Riemann invariant is extrapolated from the interior of 
the computational domain. 

For subsonic outflow the circumferential and radial 
flow angles, total pressure, and the total temperature 
are extrapolated from the interior of the computational 
domain. The total -to-static pressure ratio is specified 
at mid-span of the computational exit and the pressure 
at all other radial locations at the exit is obtained by 
integrating the equation for radial equilibrium. 
Periodicity is enforced along the outer boundaries of 
the H-grids in the circumferential direction. 

At solid surfaces the relative velocity is set to zero, 
the normal derivative of the pressure is set to zero, 
and the surfaces are assumed to be adiabatic. 


GEOMETRY AND FLOW CONDITIONS 

The two-stage supersonic turbine configuration is 
typical of those proposed for a reusable launch 
vehicle. The experimental rig is a 70% scale version 
of the actual turbine and has 12 first-stage vanes, 30 
first-stage rotors, 73 second-stage vanes and 56 
second-stage rotors (see Fig. 3). In the current effort, 
a 15-vane/30-rotor/75-vane/60-rotor airfoil 
approximation has been made, resulting in a 
simulation blade count of 1/2/5/4. To keep the pitch- 
to-chord ratio (blockage) constant, the first-stage 
vanes were scaled by a factor of 12/15, the second- 
stage vanes were scaled by a factor of 73/75 and the 
second-stage rotors were scaled by a factor of 56/60. 
The tip clearance in the first- and second-stage rotors 
was set at the design value of approximately 2% of 
the respective rotor heights. 

The total number of grid points used to discretize the 
turbine was 4,139,957. Figure 2 illustrates an x-y 
view of the grids at midspan, where every other grid 
point in each coordinate direction has been removed 
for clarity. The average value of y+, the non- 
dimensional distance of the first grid line above the 
surface was approximately 1.0 for the airfoil surfaces 
and 1.5 for the endwall surfaces. 

The simulations were run on 48 400-MHz R 12000 
processors of an SGI Origin2000. Each simulation 
was run for 15 global cycles (which constitutes on 
complete rotor revolution) at 22,000 time steps per 
cycle. A global cycle is defined as the time it takes for 
the two first-stage rotor blades to pass by one first- 
stage vane airfoil. The value of 22,000 iterations per 
cycle was chosen to resolve all the expected 
frequencies of interest. A time step required an 
average of 8.0 seconds CPU time on each of 48 (24 


MPI processes, with 2 OpenMP threads per MPI 
process) processors. The time periodicity of the 
solutions was determined by interrogating pressure 
traces at specific points along the airfoil surfaces. 

RESULTS 

Simulations have been performed at the design flow 
conditions, as well as four off-design flow conditions. 
The pressure ratios and rotational speeds of all five 
test cases are presented in Fig. 4 and tabulated in 
Tables 1 to 5. The four off-design flow conditions 
were chosen to lie on the bounds of the test envelope. 
Tables 1 to 5 also contain flow and performance 
quantities from the meanline and Navier-Stokes 
(CFD) analyses. 

Turbine performance curves predicted using the 
meanline and CFD codes are shown in Figs. 5 and 6. 

In Figs. 5 and 6 the pressure ratios refer to the total 
pressure at the inlet to the first-stage vane and the 
static pressure at the exit of the second-stage rotor. 
The shapes of the curves are similar, although the 
meanline code consistently predicts higher 
efficiencies. The differences in the efficiencies are to 
be expected because there is little experimental data 
available for anchoring the meanline and CFD 
analyses for supersonic turbines at off-design 
operating conditions. Other observations from 
comparing the CFD and meanline results in Tables 1 
to 5 include: 

1) The mass flow rate predicted by the CFD 
analysis is consistently 5% greater than that 
predicted by the meanline analysis. 

2) The CFD simulations show significantly 
more underturning of the flow in the first- 
stage rotor compared to the meanline results. 
The meanline analysis yields 
correspondingly more work in the first stage. 

3) In general, except for the off-2 flow 
conditions, the CFD analysis predicts more 
work being produced by the second stage. 
The resulting work splits between the first 
and second stages are different in the two 
analyses. 

One of the goals of the experimental portion of the 
development program is to provide a detailed data set 
at off-design operating conditions for anchoring 
codes. 

The remainder of this section will focus on only two 
of the operating conditions for brevity: the design 
flow condition and the case with a pressure ratio of 
approximately 12 (designated as off-3). This off- 
design condition was chosen because it showed the 
largest deviations from the design flow conditions. 
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Figures 7 and 8 display midspan instantaneous 
absolute Mach contours from the two cases. In both 
simulations the flow enters the vane passage at a low 
subsonic Mach number. The flow accelerates through 
the vane passage, reaching sonic (choked) conditions 
at the throat. The flow becomes supersonic in the 
diverging portion of the vane passage, reaching a peak 
number of approximately 1.8. A strong expansion 
wave system emanates from the pressure side of the 
trailing edge. The expansion waves move across the 
vane passage and interact with the suction surface of 
the adjacent vane. The vane wakes convect into the 
rotor passage where they interact with the rotor bow 
shock. This interaction creates a subsonic flow region 
near the leading edge that weakens the boundary layer 
and increases the extent of the rotor suction surface 
separated flow region. This condition exists until the 
rotor passes through the vane wake. The flow 
diffuses as it exits the rotor passage, decelerating back 
to subsonic flow conditions. The fluctuations in the 
size of the separated flow region on the suction 
surface of the rotor result in incidence variations in 
the flow entering the second-stage vane. The 
incidence variations induce periodic separation on the 
pressure surface of the vane. During the turbine 
design process the axial gap between the two stages 
was varied to minimize the incidence variations. The 
flow reaccelerates in the vane passage, reaching sonic 
conditions near the passage exit. The subsequent 
shock pattern is stronger at the off-design flow 
conditions. In both cases the flow remains subsonic 
in the second-stage rotor, although the exit Mach 
number is higher at the off-design conditions because 
of the larger pressure ratio (also see Tables I and 4). 

Figures 9 and 10 show instantaneous entropy contours 
at midspan for the two cases. Both figures illustrate 
the “wagging” of the first-stage vane wake due to the 
periodic passing of the rotor blades and the interaction 
with the rotor shock system. The variations in the 
rotor suction surface separation are highlighted in the 
rotor passage. The increased separation at the off- 
design flow conditions causes underturning of the 
flow, which is also evident by comparing the time- 
averaged rotor exit flow angles in Tables 1 and 4. 
Vortex shedding from the trailing edge of the second- 
stage vane is observed at both flow conditions, 
although it appears to be stronger at the off-design 
conditions. 

Unsteady pressure envelopes for the two operating 
conditions are shown in Figs. 11 and 12. The 
predicted solutions for the first-stage vane are similar 
at both operating conditions, which is expected since 
the vane remains choked at both pressure ratios. The 
only appreciable unsteadiness on the first-stage vane 


occurs downstream of the throat on the suction 
surface. The first-stage rotor displays significant 
unsteadiness over the first 25% of the axial chord on 
the suction surface. This unsteadiness is generated 
through a combination of the passing vane wakes, the 
vane trailing edge expansion fans, and the temporal 
variations of the rotor bow shock. The amplitude of 
the unsteadiness in this region is greater at the off- 
design conditions. The remainder of the rotor surface 
(at both operating conditions) exhibits a moderate 
amount of unsteadiness, which is generated by the 
first-stage vane wakes, interaction with the potential 
fields of the upstream and downstream vane rows and 
the intermittent separation of the rotor suction-surface 
boundary layer. The effects of the suction surface 
flow separation are evident from approximately 
X/C=0.80 to the trailing edge. The pressure envelopes 
on the second-stage vane and rotor are, in general, 
similar at both flow conditions. The exception is near 
the along the suction surface of the second-stage 
rotor. The amplitude of the unsteadiness in this region 
is slightly greater at the design flow conditions. 

Unsteady pressure traces and Fourier decompositions 
at midspan of the first-stage vane trailing edge, first- 
stage rotor leading edge and second-stage vane 
trailing edge are shown in Figs. 13-24. These 
locations were chosen for presentation because they 
exhibit the greatest amount of unsteadiness. The 
pressure variations at the trailing edge of the first- 
stage vane (see Figs. 13-16) are mainly at the rotor- 
passing frequency (5206 Hz at the design conditions 
and 5468 Hz for the off-3 conditions). The 
distribution at the off-design flow conditions also 
displays significant unsteadiness at half the rotor 
passing frequency. This additional frequency may be 
a result of increased “wagging’* of the vane wakes. A 
large amount of the unsteadiness at the leading edge 
of the first-stage rotor (see Figs. 17-20) is generated 
by the periodically passing vane wakes (2608 Hz and 
2734 Hz, respectively). The unsteadiness at the wake 
passing frequency (the large peaks in Figs. 17 and 19) 
is enhanced by the fact that these wakes periodically 
eliminate the rotor leading edge bow shock. A 
significant amount of unsteadiness is also generated at 
twice the vane wake passing frequency, and is caused 
by periodic variations in the rotor bow shock system 
and the interaction of the bow shocks with the 
adjacent rotors. The frequency spectrum for the rotor 
leading edge is much richer (in both cases) than that 
observed at the vane trailing edge. Undoubtedly the 
periodic fluctuations in the extent of the rotor suction 
surface separation and variations in the strength of the 
shock system contribute to these additional 
frequencies. The pressure traces at the trailing edge of 
the second-stage (see Figs. 21-24) vane exhibit the 
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high-frequency unsteadiness associated with vortex 
shedding. The vortex shedding frequency occurred at 
approximately 110 kHz (corresponding to a Strouhal 
number of approximately 0.22). The vortex shedding 
does not occur continuously, but appears and 
disappears as a function of the position of the first- 
stage rotor. This may be a direct consequence of the 
incidence variations discussed previously. 

Radial profiles of the circumferentially and time- 
averaged absolute Mach number are shown in Figs. 
25 and 26. The shapes of the profiles are similar in the 
first-stage for both operating conditions, which is 
expected because the first-stage vane is choked in 
both cases. The Mach number profiles are also similar 
at the exit of the second-stage vane. The shapes of the 
profiles at the exit of the second-stage rotor are 
similar, but the off-design case shows higher 
quantitative levels because the flow expands to 
achieve the larger pressure ratio. In both simulations 
the footprint of the rotor tip clearance flows can be 
seen outboard of 80% span. In both simulations there 
is more spanwise distribution of the flow in the first- 
stage rotor than in any other blade row. This 
phenomenon is caused by the rotor bow shock system 
and separated flow region, which vary spatially as 
well as temporally. 

Integrated unsteady forces on the rotors are shown in 
Figs. 27-30, respectively. The unsteady forces 
(amplitude and magnitude) on the first-stage rotor are 
similar at both operating conditions. The variations in 
the tangential force occur primarily at the passing 
frequency of the first-stage vanes. The average forces 
on the second-stage rotor are greater at the off-design 
condition, which is again consistent with the larger 
pressure ratio. The force variations in the second stage 
occur across a wide range of frequencies. 

CONCLUSIONS 

A series of meanline and unsteady three-dimensional 
Navier-Stokes simulations have been conducted to 
provide pre-test performance and unsteadiness 
predictions for a two-stage supersonic turbine. The 
performance and flow quantities predicted using the 
meanline and Navier-Stokes analyses show 
reasonable agreement over a wide range of flow 
conditions. Investigation of the unsteady flow field 
revealed a complex interaction of wakes, expansion 
fans, shock waves and separated flow regions. The 
unsteady pressure fields predicted in this study were 
used to help determine the locations of the transducers 
in the experiments. The predicted results will be 
compared to the experimental data during the next 


several months, and will be presented in a future 
paper. 
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Table 1. CFD and meanline predictions at design flow conditions 
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Table 2. CFD and meanline predictions at off-1 flow conditions 



Table 3. CFD and meanline predictions at off-2 flow conditions 
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Table 4. CFD and meanline predictions at off-3 flow conditions 
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Table 5. CFD and meanline predictions at off-4 flow conditions 
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Figure 1. Meanline code validation cases 
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Figure 2. Turbine midspan computational grid 
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Figure 5. Efficiency vs. velocity ratio 
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Figure 6. Efficiency vs. rotation speed 







Figure 8. Instantaneous absolute Mach number - off-3 



Figure 7. Instantaneous absolute Mach number - design 


12 




Figure 10. Instantaneous non-dimensional entropy - off-3 
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Figure. 11. Unsteady pressure envelopes at midspan of the first stage (design at left, off-3 at right, vanes 

top, rotors bottom) 



Figure. 12. Unsteady pressure envelopes at midspan of the second stage (design at left, off-3 at right, vanes 

top, rotors bottom) 
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Figure 13. Unsteady pressure trace - midspan of first-stage vane trailing edge - design 
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Figure 14. Fourier decomposition - midspan of first-stage vane trailing edge - design 
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Unsteady pressure trace - midspan of first-stage vane trailing edge - off-3 
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Figure 16. Fourier decomposition - midspan of first-stage vane trailing edge - off-3 
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Figure 17. Unsteady pressure trace - midspan of first-stage rotor leading edge - design 
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Figure 18. Fourier decomposition - midspan of first-stage rotor leading edge - design 
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Figure 19. Unsteady pressure trace - midspan of first-stage rotor leading edge - off-3 



Figure 20. Fourier decomposition - midspan of first-stage rotor leading edge - off-3 
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Figure 21. Unsteady pressure trace - midspan of second -stage vane trailing edge - design 



Figure 22. Fourier decomposition - midspan of second-stage vane trailing edge - design 
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Figure 23. Unsteady pressure trace - midspan of second-stage vane trailing edge - off-3 



Figure 24. Fourier decomposition - midspan of second-stage vane trailing edge - off-3 
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Figure 25. Radial profiles of circumferentially and time-averaged absolute Mach number - design. 



Figure 26. Radial profiles of circumferentially and time-averaged absolute Mach number - off-3 
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Figure 27. Unsteady integrated forces on the first-stage rotor - design 



Figure 28. Unsteady integrated forces on the first-stage rotor - off-3 
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Figure 29. Unsteady integrated forces on the second-stage rotor - design 



Figure 30. Unsteady integrated forces on the second-stage rotor - off-3 
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